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Proposed even- denominator fractional quantum Hall effect (FQHE) states suggest the possibility 
of excitations with non-Abelian braid statistics. Recent experiments on wide square quantum wells 
observe even-denominator FQHE even under electrostatic tilt. We theoretically analyze these struc- 
tures and develop a procedure to accurately test proposed quantum Hall wavefunctions. We find 
that tilted wells favor partial subband polarization to yield Abelian even-denominator states. Our 
results show that tilting quantum wells effectively engineers different interaction potentials allowing 
exploration of a wide variety of even-denominator states. 

PACS numbers: 73.43.-1,71. lO.Pm 



Experimental evidence for more than seventy odd- 
denominator FQHE states [l| demonstrates the ubiquity 
of incompressible quantum Hall electron liquids in high 
mobility semiconductor quantum wells. The Laughlin 
wavefunctions [H describe a few of these states, how- 
ever, the composite fermion (CF) theory includes Laugh- 
lin states into a comprehensive framework that captures 
the essential physics of the observed lowest Landau level 
(LLL) FQHE states I,!]. 

Observations of rare even-denominator FQHE states, 
exceptions to the odd-denominator rule, occur in multi- 
coinponent systems @ and at half filling of the second 
LLl6|, i.e., filling factor 5/2. Spin, layer, or subband de- 
grees of freedom allow combinations of single-component 
odd-denominator states to yield even-denominator frac- 
tions Q- But at half filling of the second LL theoretical 
analyses [1] suggest that the Coulomb interaction hap- 
pens to favor the formation of a single component paired 
CF state, the Pfaffian Q, a profoundly distinct state. 

Proposed FQHE states possess topologically non- 
trivial excitations. The common CF states carry excita- 
tions that obey Abelian anyonic braid statistics. Braid- 
ing two Abelian excitations changes the overall wavefunc- 
tion by a phase In contrast, excitations above the 

CF paired Pfaffian state obey non-Abelian braid statis- 
tics; a braid operates on the overall wavefunction by a 
matrix. Non-Abelian excitations have potential applica- 
tion in topological quantum information processing mak- 



ing observation of non-Abelian excitations a key goal [11 



Experimental parameters (width and density) in nar- 
row quantum well samples allow only limited tunabil- 
ity. Recent diagonalization studies |12| indicate that in a 
single-layer system the Pfaffian state is likely (unlikely) 
to be the ground state for = 5/2(1/2) within the simple 
Coulomb interaction model. By chance, including finite 
thickness of realistic quantum wells marginally favors the 
formation of the Pfaffian state in the second LL 131. A 



phases and their phase boundaries requires construction 
of tunable high mobility quantum wells designed to favor 
these fragile states by engineering the effective electron- 
electron interaction. 

Recent experiments with wide square quantum wells 
demonstrate even-denominator LLL FQHE in a reg ime 
with a surprisingly strong inter- layer tunneling IJ, [15 1 
suggesting that explanations in terms of the usual 
Abelian multicomponent states need to be revisited. Do 
these observations suggest that the new samples effec- 
tively tailor the interaction to favor a LLL Pfaffian state 
[isl llij ? The competition between the (non-Abelian) 
Pfaffian and (Abelian) multicomponent states remains a 
subtle issue [l^, [l3l- Furthermore, it is currently un- 
known whether these or similar samples will favor a sec- 
ond LL Pfaffian or entirely new FQHE states. 

In this Letter we show that recent experiments observ- 
ing even-denominator FQHE [l3 . 15 1 in wide quantum 
wells occur in a regime that can favor partial subband oc- 
cupancy and therefore offer the ability to engineer a wide 
array of even-denominator FQHE states. We develop a 
protocol that combines a modified high field local density 
approximation (LDA) of quantum well subbands with 
FQHE wavefunctions in the plane of the quantum well. 
Focusing on the experimental parameters of Ref. [isj we 
argue that evidence of even-denominator FQHE in tilted 
samples arises from partially subband polarized Abelian 
FQHE states [l3| and, in turn, demonstrate the remark- 
able ability to tune among subband balanced and imbal- 
anced FQHE states. Our protocol also applies to the sec- 
ond LL where subband engineering can be used to tune 
and explore second LL even-denominator states [l7| . 

We begin with a model of square quantum wells of wide 
width where at most two subbands are populated. We 
construct a procedure to energetically minimize the fol- 
lowing three dimensional model in a strong perpendicular 
magnetic field, B: 



route towards the exploration of putative non-Abelian 



H^" = K + Vw + V + Hb, 



(1) 



where K is the electron kinetic energy, Vvv is the quantum 
well confinement potential, V = ^i^j e^/(2e|ri — rj|) is 
the three dimensional Coulomb interaction, e is the di- 
electric constant, and 77b accounts for the energy of a 
plane of rigid positive background charges. In what fol- 
lows we consider a tilted square quantum well to model a 
specific set of experimental parameters [isj as an applica- 
tion of our procedure: V-\f^{z±) = az± + Vo9{\z±\ — w/2). 
The extra potential, az±, tilts the well. Wc choose 
Vq ~ 270 meV and a well width of w = 55 nm. 

To determine the tilt parameter a we first consider the 
B = limit. To solve for a we assume a Fermi sea in 
the plane and apply usual LDA methods . Here we 
assume that the electron spins are unpolarized and use 
the Hedin-Lundqvist exchange correlation energy [20| at 
a planar density of p = 1.72 x 10^^ 
the density imbalance (7.8 x 10^° cm 
splitting measured in Ref. ^ (A^i=° = 41 K) are re- 
produced with a = 0.3 meV/nm. 

We use this value of a to connect to the high field 
experiments of Ref. [3l and now consider the high 
field, LLL limit of H^^. The planar and perpendic- 
ular coordinates separate and we assume that in the 
xy plane the electrons are in the LLL with a quenched 
kinetic energy and basis states given by: 4>m{z/lo) = 

where z 



cm . We find that 
^) and subband 



(z/lo)"' exp (-|z|V 4;g)/(/oV2^2'»m! 
and lo = yJhc/eB is the magnetic length. We also as- 
sume that the real spins are fully polarized (the experi- 
ments at total filling = 1/2 in Ref. [l^ are performed 
for B « 14T). We then use LDA with a spin polarized 
exchange correlation energy [2l| to model the electron 
wavefunction perpendicular to the plane. 

In wells of wide width the subband polarization, 7 — 
{Nq — Ni)/N^ is a key unknown connecting the energetics 
of the plane to the energetics arising from perpendicular 
coordinates. Here Nq (Ni) denotes the number of elec- 
trons in the lowest (first) subband with = Ao + A^i. 
At zero field 7 is determined by the energetics of the xy- 
plane Fermi surface [l^ but this is not obviously accurate 
in wide quantum wells in the LLL. 

To determine the ground state we minimize the to- 
tal energy per particle, i^*"'*^' = + E^^ , where E^ = 
{NoEo+NiEi)/N, the weighted sum of both subband en- 
ergies, arises from the quantum well confinement and EH^^ 
is the correlation energy due to the LLL planar Coulomb 
interaction. We calculate E*^°^^^ with the following pro- 
cedure: For each value of 7 we use LDA to compute 
the confinement energy, E;!^ — A^/2 — 7Aq]^/2, rewritten 
here in terms of the 7 dependent subband energy differ- 
ence, Agj, and an effective chemical potential, A^. The 
LDA also yields ^a, where ^0 and ^1 denote the lowest 
and first subband wavefunctions, respectively. The LDA 
output is used to construct an effective 2D model of spin 
polarized LLL fermions, = H:^ + H^^, one at each 
7. Finally, we compute E^'^ with two-component varia- 
tional FQHE wavefunctions in the subband basis at fixed 
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Figure 1: Top: Confinement energy plotted as function of 
subband polarization for parameters corresponding to tilted 
wide well samples of Ref. [1^ . The right axis shows Coulomb 
units in terms of the fixed magnetic length, I, obtained for 
u = 1/2 and p = 1.72 x 10"cm~^. Bottom: The same but for 
the subband energy diflerence (dashed) and the chemical po- 
tential offset (solid). The dotted line plots the B = subband 
energy difference. Inset: Schematic showing partial subband 
occupancy captured by an effective chemical potential. 



7 thus allowing global minimization of E^°^^^ in the space 
of competitive FQHE wavefunctions. Our procedure re- 
sults in a set of models, H°^, that are very sensitive to the 
quantum well parameters showing that experiments can 
tune through a large set of multicomponent interactions. 

To model the perpendicular degrees of freedom wc in- 
troduce an effective model that captures the subband 
energetics: 



E 



A1 

,,0) TT + nmfi) 



where ' 



4 



and cl^^ creates an electron in the 
state £,a(j)m- H'^ parameterizes interacting electrons at 
fixed 7 in the perpendicular direction as non-interacting 
fermions, with an effective chemical potential. The first 
term is the subband splitting commonly used in models 
of very wide wells but the second term is added to capture 
the energetics of partial subband occupation in wide wells. 
The inset of Fig. ([1]) shows that A^ > favors partial 
subband occupancy. We obtain the very wide well limit 
(two fully-occupied subbands) for A^ ^ Aqi and the 
narrow well limit (one subband occupied) with A^ ^ 
Aoi. 

In the top panel of Fig. ([T]) we see that the lowest 
energy is obtained for subband polarization 7 « 0.59. 
The large upturn near 7 = arises because the well 
tilt strongly penalizes occupancy of the second subband. 
The upturn near 7 = 1 arises from the Coulomb cost 
of putting all charges in the lowest subband. The bot- 
tom panel of Fig. ([T]) plots A^^ and Aqj^ versus 7 to show 
that they vary in comparison to the subband splitting 



3 



-0.2 



-0.25 




-4 -2 2 4 
z^[nm] / 



Mg=2.55 




(3/7,1/511) -pO.S 

--- (1/2,1/212) it0.5 
. -. (l/2,l/2l2)pjit0.5 



Figure 2: In-plane correlation energies of several states, 
Eqs. ([3]) and Q, at 7 = 0.5 plotted versus the inter-subband 
interaction parameter R. The arrow indicates the parame- 
ter relevant for the tilted well samples of Ref. [3| where the 
(3/7, 1/5|1) state is shown to have the lowest correlation en- 
ergy at 7 = 0.5. Inset: Normalized lowest (solid) and first 
(dashed) subband wavefunctions plotted versus position in 
the tilted quantum well for several 7 between and 1. In the 
limit 7 — >■ the subbands localize to mimic layer states. 



obtained from the zero field calculation, A^j^-°. The 
non-linearity demonstrates a Stoner-like dependence on 
the subband-pseudospin occupancy that arises from the 
competition between the kinetic energy, well tilt, and 
Coulomb interaction along the direction perpendicular 
to the plane. 

We now use output from our 7 dependent LDA cal- 
culation, ^cr, to construct an effective interaction within 
the plane. The LLL interaction contains matrix elements 
of the form: VZ\''^^p^^ = (mi, mal^M (r)|m3, ^4). 
The effective planar interaction is given by: Vs^^y (r) = 
(C<Ti : ^(72 1 y I ^0-3 , (,ai ) ■ The tilt breaks the symmetry in the 
interaction, e.g., Vioio 7^ Vbioi- The inset of Fig. [5] plots 
the subband wavefunctions. The tilt localizes the sub- 
band states on either side of the well. 

For the tilted wide well we find that we can ignore 
Viioo and Vooii in the regime < 7 < 0.7 and 7 = 1 to 
yield a simple effective model: 



-'-'7 ~ -'-'7 ^ 9 {m} ^mi(T^m2cr'^m4cr'^m3CT- ) 

The off-diagonal terms give Haldane pseudopotentials 
(2^ that are smaller by more than a factor of 7 than the 
diagonal terms. We assume that diagonal repulsion dom- 
inates the inter-subband energetics and that the omission 
of these small off-diagonal terms docs not change the en- 
ergetic ordering of trail wavefunctions computed below 



using Eq. ([2]). 

We model the planar interaction computed with the 
subband wavefunctions by a modified Zhang-Das Sarma 
23| potential: 



V 



MZDS 



where the parameters d, c, and W are determined by fit- 
ting lowest pseudopotentials of y^^^DS those of the 
interaction computed using the subband wavefunctions. 
Wc have found that y'^ZDS gjvQg good fits for a variety of 
quantum wells in the subband basis. For non-tilted wide 
wells the c term must be retained. We find largest uncer- 
tainty (< 5%) in the m = 1 pseudopotential. We have 
checked, by varying fitting parameters, that 5% varia- 
tion in pseudopotentials does not qualitatively alter our 
results. 

For the tilted well sample of Ref. [lij we focus on 
the regime < 7 < 0.7 and 7 = 1 and keep only the 
Vaa'aa' tcrms. Wc then require only three fitting param- 
eters: do(7) = c^oooo, ^1(7) = chilli, and R{j) = dioio 
to approximately match all pseudopotentials. For ex- 
ample, we find the lowest E^°^^^ to occur at: do (0-5) = 
1.16^0,^1(0.5) = 1.37^0, and i?(0.5) = 2.55^0- 

Equation ^ resembles models of bilaycrs and very 
wide wells but there arc crucial differences: (i) There 
is a term that favors partial subband polarization, 
(ii) All tcrms arc functions of the subband polarization, 
7, (iii) The interaction within each subband differs, c.g, 
^1111 7^ Voooo, and (iv) In the absence of tilt we must 
keep off-diagonal terms of the form V^aa'a'- 

We now construct competitive two-component varia- 
tional wavefunctions expected to minimize the energy 
of Eq. ^ and therefore H^^ . The single component 
LLL CF wavefunctions at filliiig n/(2pn± 1), V'n/(2pn±i)j 
are given in the literature |j| . The single component 
CF wavefunctions can be generalized to capture two- 
component states [1^ : 



{i'o,i'i\m} 



= U^^J - Wrr^hoiUkW^AM], (3) 



where the fully antisymmetric wave function is a sin- 
gle component state at filling factor V. Halperin's wave 
functions are obtained as special cases for Vq = l/m' 
and Vi = l/m" . In order to ensure that the electrons 
of each component occupy the same area, TVq and Ni 
must be related by Nov^"^ -t- mA'^i — iViFJ^^ + mNo, thus 



Some of 



the 7 = states of Eq. ^ have been shown to be energet- 
ically competitive with favorable ground states at several 
different filling factors in bilayer systems (small tunnel- 
ing and equal subband population) [l3|. Table H] lists a 
larger set. The constituent wavefunctions, ipv; with the 
lowest energy gap yield an upper-bound for the energy 
gap of each multicomponcnt state. For example, the gap 
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Table I: Several possible states, Eq. Q, at half filling. The 
last column, polarization, shows that two of the incompress- 
ible states are locked at partially polarized configurations. 



of the (3/7, 1/5|1) state is no larger than that of the -01/5 
state while (1/2, 1/2|2) is a gapless CF Fermi sea We 
also include the multicomponent Pfafhan states in our 
comparison: 



n( 



-u;,)'"7/.pf[{zfc}]V'Pf[{u;J], (4) 



where ippf is the Pfaffian wavefunction Q. 

We compute the ground state energies of several can- 
didate wavef unctions in the subband basis at v = 1/2 
for the parameters of Ref. [l^ using variational Monte 
Carlo. We use the spherical geometry Q to compute en- 
ergies in finite size systems and extrapolate our results to 
the thermodynamic limit by adding the background en- 
ergy: E^, =' - [N?F{di) + N^Fjdo) + 2NoNiF {R)] /2, 
where F{x) = {e^ /2elo){-x/ + and Rs 

is the radius of the sphere in units of the magnetic length. 

Fig. [2] shows a representative comparison of E^^ for 
three of the lowest energy states at 7 = 0.5. Here the 
parameter R is varied to test the robustness of the lowest 
energy state. The Monte Carlo error bars are smaller 
than the line width. The lowest energy state at 7 = 0.5 
is found to be the incompressible (3/7, 1/5|1) state. The 
compressible Fermi sea state, (1/2,1/2|2), is nearby in 
energy. The Pfaffian state becomes lower in energy for 
the unphysical regime of R < lo- 

We have compared wavefunctions at several subband 
polarizations to globally minimize the total energy. We 
expect partially subband polarized states to be favored 
by the competition between the Coulomb interaction 
along the direction perpendicular to the plane and the 
quantum well tilt. We find that the small gain in E-^ fa- 
vors the (3/7, 1/5|1) state over all others in a parameter 
window corresponding to the parameters relevant for the 
tilted well experiments of Ref. 15|, a = 0.3 meV/nm, see 
Table [ni We also conclude that small changes in a will 
favor the partially subband polarized compressible state, 
(1/2,1/2|2). 

Our procedure provides the following physical picture 
of the experiments in Ref. [3l at = 1/2. A symmetric 
wide well favors the compressible (1/2, 1/2 1 2) state with 
7 « (or the incompressible 7 = (1/3,1/3|1) state 
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Table II: Comparison of planar correlation energies (third col- 
umn) and confinement energies (fourth column) for several 
V = 1/2 states. The energies were computed for parameters 
relevant for the tilted well sample of Ref. [l^ . The lowest to- 
tal energy is found for the incompressible (3/7, 1/5|1) state. 



depending on p and w). As the density is made more 
asymmetric via a tilt the CF Fermi sea state continu- 
ously depopulates its higher subband (it is a subband- 
paramagnet). Near a « 0.3 meV/nm each subband be- 
comes localized on either side of the well to approximate 
layer-like behavior but with asymmetric parameters in 
the subband basis. The incompressible asymmetric state, 
(3/7, 1/5|1), is then favored in a narrow parameter win- 
dow. We find that this state is responsible for the indi- 
cations of FQHE observed under tilt in Ref. [l3| • As the 
state is made more asymmetric the higher subband state 
becomes further depopulated to again favor (1/2, 1/2|2). 
Our analysis shows that the samples of Ref. [l5[ allow 
tuning among partially subband polarized multicompo- 
nent states. 

We have shown that quantum wells with two active 
subbands can exhibit a wide variety of partially subband 
polarized, even-denominator FQHE states. Our proce- 
dure can be used to find parameter regimes where incom- 
pressible partially subband polarized FQHE states arise 
as a function of sample density, width, and tilt (e.g., a 



mixed paired state [24 
for recently observed 



]. E^ Q also provides candidates 
LLL = 1/4 features (e.g., 
(3/13, 1/7|3)), and the second LL FQHE where tilt can 
be used to tune between a variety of subband polariza- 
tions thereby allowing one to engineer a larger class of 
interaction potentials. Furthermore, Eq. ([3|) offers sev- 
eral candidate partially polarized real-spin FQHE states. 
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